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Abstract 

<D 

We present hadron-level predictions from the Monte Carlo generator Cascade and nu- 
■ merical level calculations of beauty quark and inclusive 6-jet production in the framework 
of the /c^-factorization QCD approach for CERN LHC energies. The unintegrated gluon 
densities in a proton are determined using the CCFM evolution equation and the Kimber- 
Martin-Ryskin (KMR) prescription. We study the theoretical uncertainties of our calcula- 
tions and investigate the effects coming from parton showers in initial and final states. Our 
predictions are compared with the recent data taken by the CMS collaboration. 
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1 Introduction 

Beauty production at high energies is subject of intense studies from both theoretical 
and experimental points of view since events containing b quarks present an important 
background to many of the searches at the LHC. From the theoretical point, the dominant 
production mechanism is believed to be quark pair production through the gluon-gluon 
fusion subprocess and therefore these processes provide an opportunity to test the different 
predictions based on Quantum Chromodynamics (QCD). The present note is motivated by 
the recent measurements [H[2] b production performed by the CMS collaboration. The b- 
quark cross sections have been presented [1] as a function of the muon transverse momentum 
and pseudorapidity at y/s = 7 TeV. It was observed that the data tend to be higher than 
the MC@NLO [3111] predictions. On the other hand the measurements of inclusive 6-jet 
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production cross sections [2] are reasonably well described by MC@NLO. In addition to the 
comparison of CASCADE with data in |1] we present here futher studies. 

In the framework of the /c^-factorization approach of QCD [5], which is of primary con- 
sideration in this note, a study of the heavy quark production has been done (for previous 
results see |6HT2]). In our previous study (12] we show a good agreement between the Teva- 
tron data on the b quarks, bb di-jets, B + and several D mesons (or rather muons from their 
semileptonic decays) production with the predictions coming from fc^-factorization and we 
investigated the role of initial and final state parton showers. Based on these results, we 
give here a systematic analysis of the recent CMS measurements [TJ[2] in the framework of 
kx~ factorization. As done in [12], we produce the calculations in two ways: we will per- 
form numerical parton-level calculations (labeled as LZ) as well as calculations with the full 
hadron level Monte Carlo event generator Cascade [13J and compare both with the mea- 
sured cross sections of heavy quark production. In this way we will investigate the influence 
of parton showers in initial and final states for the description of the data. Additionally 
we study different sources of theoretical uncertainties, i.e. uncertainties connected with the 
gluon evolution scheme, heavy quark mass, hard scale of partonic subprocess and the heavy 
quark fragmentation functions. 

The outline of our paper is the following. In Section 2 we recall very shortly the basic 
formulas of the /^-factorization approach with a brief review of calculation steps. In Section 3 
we present the numerical results of our calculations and a discussion. Section 4 contains our 
conclusions. 

2 Theoretical framework 

In the present analysis we follow the approach described in the earlier publication |12j . 
For the reader's convenience, we only briefly recall here main points of the theoretical scheme. 

The cross section of heavy quark hadroproduction at high energies in the /^-factorization 
approach is calculated as a convolution of the off-shell (i.e. /^-dependent ) partonic cross 
section a and the unintegrated gluon distributions in a proton. It can be presented in the 
following form: 



where A(x,k^, fi 2 ) is the unintegrated gluon distribution in a proton, \M.(g*g* — > QQ)\ 2 is 
the off-shell (i.e. depending on the initial gluon virtualities k 2 T and k| T ) matrix element 
squared and averaged over initial gluon polarizations and colors, and s is the total center-of- 
mass energy. The produced heavy quark Q and anti-quark Q have the transverse momenta 
Pit and p2T and the center-of-mass rapidities y± and y<i- The initial off-shell gluons have a 
fraction x\ and x^ of the parent protons longitudinal momenta, non-zero transverse momenta 
kir and k 2 T (k^ T = — k\ T ^ 0, k 2 T = —k\ T ^ 0) and azimuthal angles 0i and 2 - The 
analytic expression for the \Ai(g*g* — > QQ)\ 2 can be found, for example, in [5|,[9]. 

The unintegrated gluon distribution in a proton A(x, k^, fi 2 ) in (1) can be obtained from 
the analytical or numerical solution of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) jH] or 
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Ciafaloni-Catani-Fiorani-Marchesini (CCFM) [15] equations. As in [12], in the numerical cal- 
culations we have tested a few different sets, namely CCFM AO (BO) [16J and KMR [T7] ones. 
The input parameters in both CCFM-evolved gluon densites have been fitted [16] to describe 
the proton structure function F 2 (x, Q 2 ). The difference between AO and BO sets is connected 
with the different values of soft cut and width of the intrinsic distribution. A reasonable 
description of the F 2 data can be achieved [16] by both these sets. To evaluate the unin- 
tegrated gluon densities in a proton A(x, k^, /i 2 ) we apply also the Kimber-Martin-Ryskin 
(KMR) approach [17]. The KMR approach is a formalism to construct the unintegrated 
parton (quark and gluon) distributions from the known conventional parton distributions. 
For the input, we have used the standard GRV 94 (LO) [18] (in LZ calculations) and MRST 
99 [20] (in Cascade) sets. 

3 Numerical results 

The unintegrated gluon distributions to be used in the cross section (1) depend on the 
renormalization and factorization scales \ir and \ip. Following [12], in the numerical calcula- 
tions we set fi 2 R = m 2 Q + (p 2 T + p 2 T )/2, n 2 F = s + Q 2 -., where Qt is the transverse momentum 
of the initial off-shell gluon pair, m c = 1.4 ± 0.1 GeV, = 4.75 ± 0.25 GeV. We use the 
LO formula for the coupling a s (fi 2 R ) with rif = 4 active quark flavors at Aqcd = 200 MeV, 
such that a s (Ml) = 0.1232. 

We begin the discussion by presenting our results for the muons originating from the 
semileptonic decays of the b quarks. The CMS collaboration has measured [T] the transverse 
momentum and pseudorapidity distributions of muons from 6-decays. The measurements 
have been performed in the kinematic range p 1 ^ > 6 GeV and \r]^\ < 2.1 at the total center- 
of-mass energy y'i = 7 TeV. To produce muons from 6-quarks in the LZ calculations, we first 
convert 6-quarks into B mesons using the Peterson fragmentation function with default value 
eb = 0.006 and then simulate their semileptonic decay according to the standard electroweak 
theory. The branching of b — > \l as well as the cascade decay b — > c — > \i are taken into 
account with the relevant branching fractions taken from [22J. The predictions of the LZ and 
Cascade calculations are shown in Figs. 1 and 2 in comparison with the CMS data. We 
find a good description of the data when using the CCFM-evolved (AO) gluon distribution in 
LZ calculations although the Cascade curves tend to lie slightly below the data at central 
rapidities. The predictions between the LZ and Cascade calculations agree well at parton 
level. The observed difference between them in Figs. 1 and 2 is due to missing parton shower 
effects in the LZ calculations. The influence of such effects is demonstrated in Fig. 3, where 
we show separately the results of our Cascade calculations without parton shower, with 
only initial state, with only final state and with both initial and final state parton showers. 
One can see that without initial and final state parton showers, the Cascade predictions 
are very close to the LZ ones. The similar situation was pointed out in [12] at the Tevatron 
energies. 

To investigate the dependence of our predictions on the quark-to-hadron fragmentation 
function, we repeated our calculations with the shifted value of the Peterson shape parameter 
€b = 0.003, which is is often used in the NLO pQCD calculations. Additionally, we have 
applied the non-perturbative fragmentation functions which have been proposed in [23TI25] . 
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CMS data [/ib] 


1.32 ±0.01 (stat) ±0.30 (syst) ±0.15 (lumi) 


AO (LZ/Cascade) 


1.31/0.96 


BO (LZ/Cascade) 


0.98/0.72 


KMR (LZ/Cascade) 


0.91/0.59 


MC@NLO [I] 


0.95 


Pythia pQ 


1.9 



Table 1: The inclusive b-quark production cross section in pp collisions at y'i = 7 TeV. 

The input parameters were determined [2UI25] by a fit to LEP data. The results of our 
calculations are shown in Fig. 4. For illustration, we used here the CCFM AO gluon density. 
We find that the predicted cross sections in the considered kinematic region are larger for 
smaller values of the parameter or if the fragmentation function from [23-25J is used. 
Thus, the CMS data points lie within the band of theoretical uncertainties. The results 
obtained here (see Fig. 1) with the CCFM B0 and KMR gluon densities (but also with AO 
density as shown in the CMS paper) are rather close to the MC@NLO ones (not shown) and 
underestimate the data by a factor of 1.6. 

The visible cross sections of decay muons from 6-decays are listed in Table 1 in comparison 
with the CMS data p]. In Table 2 the systematic uncertainties of our calculations are 
summarized. To estimate the uncertainty coming from the renormalization scale fiR, we 
used the CCFM set A0± and AO— instead of the default density function AO. These two 
sets represent a variation of the scale used in a s in the off-shell matrix element. The A0± 
stands for a variation of 2 fi R , while set AO— reflects [Ir/I. We observe a deviation of 
roughly 13% for set A0±. The uncertainty coming from set AO— is generally smaller and 
negative. The dependence on the 6-quark mass is investigated by variation of 6-quark mass 
of rrib = 4.75 GeV by ±0.25 GeV. The calculated 6-quark cross sections vary by ~ ±6%. 

The CMS collaboration has measured [2] the double differential cross sections da/dydpr 
of inclusive 6-jet production at the yfs = 7 TeV. The measurements have been determined in 
four 6-jet rapidity regions, namely \y\ < 0.5, 0.5 < \y\ < 1, 1 < \y\ < 1.5 and 1.5 < \y\ < 2. 
Our predictions are shown in Figs. 5 and 6 and compared to the CMS data. In the Cascade 
calculations the 6-jets are reconstructed with the anti-fc t cone algorithm [26J (using the 
Fastjet package [27] 128]) with radius R > 0.5. In contrast with the decay muon cross 
sections, the predictions based on the CCFM and KMR gluons are very similar to each 
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Source 


a(pp — > 6 + A — )• /i + A', pj, > 6 GeV, \if\ < 2.1) 


CCFM set AO 


0.96 fib 


CCFM set A0+ 


+ 13% 


CCFM set A0- 


-2% 


m b = 5.0 GeV 


-7% 


m h = 4.5 GeV 


+6% 


, _ n nno 
Cfr — u.uuo 


+9% 


Total 


+ 17% 
=■=7% 



Table 2: Systematic uncertainties for beauty total cross section in pp collisions at y/s = 7 TeV 
obtained with Cascade. 

other. The reasonable description of the data is obtained by all unintegrated gluon densities 
under consideration. 

Finally, we would like to point out the role of non-zero gluon transverse momentum kx 
in the off-shell matrix elements (see Figs. 7 and 8). In these Figs, the solid histograms cor- 
respond to the results obtained according to the master formula (1). The dotted histograms 
are obtained by using the same formula but without virtualities of the incoming gluons in 
partonic amplitude and with the additional requirement k^ 2T < n R . We find that the non- 
zero gluon transverse momentum in the hard matrix element is important for the description 
of data. The similar situation was pointed out in [12] at the Tevatron energies. It means, 
that the high region is important, and only when including the highk-r tail the results 
are similar to NLO predictions. 

4 Conclusions 

In this note we analyzed the first data on the beauty production in pp collisions at LHC 
taken by the CMS collaboration. Our consideration is based on the /c;r-factorization approach 
supplemented with the CCFM-evolved unintegrated gluon densities in a proton. The analysis 
covers the total and differential cross sections of muons originating from the semileptonic 
decays of beauty quarks as well as the double differential cross sections of inclusive 6-jet 
production. Using the full hadron-level Monte Carlo generator Cascade, we investigated 
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the effects coming from the parton showers in initial and final states. Different sources of 
theoretical uncertainties have been studied. 

Our LZ predictions with the default set of parameters agree with the data. The Cascade 
predictions tend to slightly underestimate the data at central rapidities but the data points 
still lie within the band of theoretical uncertainties. In this case the overall description of the 
data at a similar level of agreement as in the framework of NLO collinear QCD factorization. 
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Figure 1: The pseudorapidity distributions of muons arising from the semileptonic decays 
of beauty quarks. The first column shows the LZ numerical results while the second one 
depicts the Cascade predictions. The solid, dashed and dotted histograms correspond to 
the results obtained with the CCFM AO, B0 and KMR unintegrated gluon densities. The 
kinematic cuts applied are described in the text. The experimental data are from CMS [I]. 
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Figure 2: The transverse momentum distributions of muons arising from the semileptonic 
decays of beauty quarks. The first column shows the LZ numerical results while the second 
one depicts the Cascade predictions. Notation of all histograms is the same as in Fig. 1. 
The kinematic cuts applied are described in the text. The experimental data are from 
CMS PQ. 
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Figure 3: Parton shower effects in the pseudo-rapidity and transverse momentum distribu- 
tions of the muons. The four lines represent full parton shower (solid line), no parton shower 
(dashed line), initial state parton shower (dashed dotted line) and final state parton shower 
(dotted line). The experimental data are from CMS pp. 
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Figure 4: The dependence of our predictions on the fragmentation scheme. The solid, dashed 
and dash-dotted histograms correspond to the results obtained using the Peterson fragmen- 
tation function with e& = 0.006, % = 0.003 and the non-perturbative fragmentation functions 
from [23TE5] . respectively. We use CCFM-evolved (AO) gluon density for illustration. The 
experimental data are from CMS [lj. 
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Figure 5: The double differential cross sections da/dydpx of inclusive 6-jet production as a 
function of pt in different y regions calculated at y/s = 7 TeV (LZ predictions). Notation of 
all histograms is the same as in Fig. 1. The experimental data are from CMS [2]. 
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Figure 6: The double differential cross sections da/dydpr of inclusive 6-jet production as 
a function of pt in different y regions calculated at yfs = 7 TeV (Cascade predictions). 
Notation of all histograms is the same as in Fig. 1. The experimental data are from CMS [2]. 
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Figure 7: Importance of non-zero transverse momentum of incoming gluons in open b quark 
production at the LHC The solid histograms correspond to the results obtained according 
to the master formula (1). The dotted histograms are obtained by using the same formula 
but now we switch off the virtualities of both incoming gluons in partonic amplitude and 
apply an additional requirement k^ 2T < fi 2 R . We have used here the CCFM AO gluon for 
illustration. The experimental data are from CMS pp. 
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Figure 8: Importance of non-zero transverse momentum of incoming gluons in 6-jet produc- 
tion at the LHC Notation of all histograms is the same as in Fig. 7. The experimental data 
are from CMS 0. 
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